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ABSTRACT 

WMAP precision data enables accurate testing of cosmological models. We find that 
the emerging standard model of cosmology, a flat A— dominated universe seeded by a 
nearly scale-invariant adiabatic Gaussian fluctuations, fits the WMAP data. For the 
WMAP data only, the best fit parameters are h = 0.72 ± 0.05, Q b h 2 = 0.024 ± 0.001, 
n m h 2 = 0.14 ± 0.02, r = 0.166l[j;$?, n s = 0.99 ± 0.04, and a 8 = 0.9 ± 0.1. With 
parameters fixed only by WMAP data, we can fit finer scale CMB measurements and 
measurements of large scale structure (galaxy surveys and the Lyman a forest). This 
simple model is also consistent with a host of other astronomical measurements: its 
inferred age of the universe is consistent with stellar ages, the baryon/photon ratio is 
consistent with measurements of the [D]/[H] ratio, and the inferred Hubble constant 
is consistent with local observations of the expansion rate. We then fit the model 
parameters to a combination of WMAP data with other finer scale CMB experiments 
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(ACBAR and CBI), 2dFGRS measurements and Lyman a forest data to find the model's 
best fit cosmological parameters: h = OJlig;^, £l b h 2 = 0.0224 ± 0.0009, ft m /i 2 = 
0.135±g;gg|, r = 0.17 ± 0.06, n s (0.05 Mpc" 1 ) = 0.93 ± 0.03, and cr 8 = 0.84 ± 0.04. 
WMAP's best determination of r = 0.17 ± 0.04 arises directly from the TE data and 
not from this model fit, but they are consistent. These parameters imply that the age 
of the universe is 13.7 ± 0.2 Gyr. With the Lyman a forest data, the model favors but 
does not require a slowly varying spectral index. The significance of this running index 
is sensitive to the uncertainties in the Lyman a forest. 

By combining WMAP data with other astronomical data, we constrain the geometry 
of the universe: Qtot = 1-02 ± 0.02, and the equation of state of the dark energy, 
w < —0.78 (95% confidence limit assuming w > — 1.). The combination of WMAP and 
2dFGRS data constrains the energy density in stable neutrinos: Vt v h? < 0.0076 (95% 
confidence limit). For 3 degenerate neutrino species, this limit implies that their mass 
is less than 0.23 eV (95% confidence limit). The WMAP detection of early reionization 
rules out warm dark matter. 

Subject headings: cosmic microwave background — cosmology: observations — early 
universe 

1. INTRODUCTION 

Over the past century, a standard cosmological model has emerged: With relatively few pa- 
rameters, the model describes the evolution of the Universe and astronomical observations on scales 
ranging from a few to thousands of Megaparsecs. In this model the Universe is spatially flat, homo- 
geneous and isotropic on large scales, composed of radiation, ordinary matter (electrons, protons, 
neutrons and neutrinos), non-baryonic cold dark matter, and dark energy. Galaxies and large-scale 
structure grew gravitationally from tiny, nearly scale-invariant adiabatic Gaussian fluctuations. 
The Wilkinson Microwave Anisotropy Probe (WMAP ) data offer a demanding quantitative test 
of this model. 

The WMAP data are powerful because they result from a mission that was carefully designed 
to limit systematic measurement errors (Bennett et al. 2003a, b; Hinshaw et al. 2003b). A critical 
element of this design includes differential measurements of the full sky with a complex sky scan 
pattern. The nearly uncorrelated noise between pairs of pixels, the accurate in-flight determination 
of the beam patterns (Page et al. 2003c, a; Barnes et al. 2003), and the well-understood properties 
of the radiometers (Jarosik et al. 2003a,b) are invaluable for this analysis. 

Our basic approach in this analysis is to begin by identifying the simplest model that fits the 
WMAP data and determine the best fit parameters for this model using WMAP data only without 
the use of any significant priors on parameter values. We then compare the predictions of this 
model to other data sets and find that the model is basically consistent with these data sets. We 
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then fit to combinations of the WMAP data and other astronomical data sets and find the best fit 
global model. Finally, we place constraints on alternatives to this model. 

We begin by outlining our methodology (§2). Verde et al. (2003) describes the details of the 
approach used here to compare theoretical predictions of cosmological models to data. In §3, we 
fit a simple, six parameter ACDM model to the WMAP data-set (temperature-temperature and 
temperature-polarization angular power spectra). In §4 we show that this simple model provides 
an acceptable fit not only to the WMAP data, but also to a host of astronomical data. We use the 
comparison with these other datasets to test the validity of the model rather than further constrain 
the model parameters. In §5, we include large scale structure data from the 2dF Galaxy Redshift 
Survey (2dFGRS, Colless et al. (2001)) and Lyman a forest data to perform a joint likelihood 
analysis for the cosmological parameters. We find that the data favors a slowly varying spectral 
index. This seven parameter model is our best fit to the full data set. In §6, we relax some of the 
minimal assumptions of the model by adding extra parameters to the model. We examine non-flat 
models, dark energy models in which the properties of the dark energy are parameterized by an 
effective equation of state, and models with gravity waves. By adding extra parameters we introduce 
degenerate sets of models consistent with the WMAP data alone. We lift these degeneracies by 
including additional microwave background data-sets (CBI, ACBAR) and observations of large- 
scale structure. We use these combined data sets to place strong limits on the geometry of the 
universe, the neutrino mass, the energy density in gravity waves, and the properties of the dark 
energy. In §7, we note an intriguing discrepancy between the standard model and the WMAP data 
on the largest angular scales and speculate on its origin. In §8, we conclude and present parameters 
for our best fit model. 



2. B AYES IAN ANALYSIS OF COSMOLOGICAL DATA 

The basic approach of this paper is to find the simplest model consistent with cosmological 
data. We begin by fitting a simple six parameter model first to the WMAP data and then to 
other cosmological data sets. We then consider more complex cosmological models and evaluate 
whether they are a better description of the cosmological data. Since Komatsu et al. (2003) found 
no evidence for non-Gaussianity in the WMAP data, we assume the primordial fluctuations are 
Gaussian random phase throughout this paper. For each model studied in the paper, we use a 
Monte Carlo Markov Chain to explore the likelihood surface. We assume flat priors in our basic 
parameters, impose positivity constraints on the matter and baryon density (these limits lie at such 
low likelihood that they are unimportant for the models. We assume a flat prior in r, the optical 
depth, but bound r < 0.3. This prior has little effect on the fits but keeps the Markov Chain out of 
unphysical regions of parameter space. For each model, we determine the best fit parameters from 
the peak of the N-dimensional likelihood surface. For each parameter in the model we also compute 
its one dimensional likelihood function by marginalizing over all other parameters; we then quote 
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the (1-dimensional) expectation value 12 as our best estimate for the parameter: 

(cti) = J d N aC(a)ai, (1) 

where a denotes a point in the N-dimensional parameter space (in our application these are points 
-sets of cosmological parameters- in the output of the Markov Chain), C denotes the likelihood 
(in our application the "weight" given by the chain to each point). The WMAP temperature (TT) 
angular power spectrum and the WMAP temperature-polarization (TE) angular power spectrum 
are our core data sets for the likelihood analysis. Hinshaw et al. (2003b) and Kogut et al. (2003) 
describe how to obtain the temperature and temperature-polarization angular power spectra re- 
spectively from the maps. Verde et al. (2003) describes our basic methodology for evaluating the 
likelihood functions using a Monte Carlo Markov Chain algorithm and for including data-sets other 
than WMAP in our analysis. In addition to WMAP data we use recent results from the CBI 
(Pearson et al. 2002) and ACBAR (Kuo et al. 2002) experiments. We also use the 2dFGRS mea- 
surements of the power spectrum (Percival et al. 2001) and the bias parameter (Verde et al. 2002), 
measurements of the Lyman a power spectrum (Croft et al. 2002; Gnedin & Hamilton 2002), su- 
pernova la measurements of the angular diameter distance relation (Garnavich et al. 1998; Riess 
et al. 2001), and the Hubble Space Telescope Key Project measurements of the local expansion rate 
of the universe (Freedman et al. 2001). 



3. POWER LAW ACDM MODEL AND THE WMAP DATA 

We begin by considering a basic cosmological model: a flat Universe with radiation, baryons, 
cold dark matter and cosmological constant, and a power-law power spectrum of adiabatic primor- 
dial fluctuations. As we will see, this model does a remarkably good job of describing WMAP TT 
and TE power spectra with only six parameters: the Hubble constant h (in units of 100 km/s/Mpc), 
the physical matter and baryon densities w m = 0, m h 2 and Wb = ftbh 2 , the optical depth to the 
decoupling surface, r, the scalar spectral index n s and A, the normalization parameter in the 
CMBFAST code version 4.1 with option UNNORM. Verde et al. (2003) discusses the relationship 
between A and the amplitude of curvature fluctuations at horizon crossing, |Ai?| 2 = 2.95 x 10~ 9 A 
In §4, we show that this model is also in acceptable agreement with a wide range of astronomical 
data. 

This simple model provides an acceptable fit to both the WMAP TT and TE data (see Figure 
1 and 2). The reduced 13 Xeff f° r ^ un ^ is 1.066 for 1342 degrees of freedom, which has a 



In a Monte Carlo Markov Chain, it is a more robust quantity than the mode of the a posteriori marginalized 
distribution. 

13 Here, xtff = — 2 In C and v is number of data minus the number of parameters. We have used 100,000 Monte Carlo 
realization of the WMAP data with our mask, noise and angle-averaged beams and found that the (— 21n£/^) = 1 
for the simulated temperature data. 
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probability of ~ 5%. For the TT data alone, Xeff/ V = 1-09, which for 893 degrees of freedom 
has a probability of 3%. Most of the excess xlff * s due to the inability of the model to fit sharp 
features in the power spectrum near / ~ 120, the first TT peak and at I ~ 350. In Figure 4 we 
show the contribution to xtff P er rnultipole. The overall excess variance is likely due to our not 
including several effects, each contributing roughly 0.5 — 1% to our power spectrum covariance near 
the first peak and trough: gravitational lensing of the CMB (Hu 2001), the spatial variations in 
the effective beam of the WMAP experiment due to variations in our scan orientation between the 
ecliptic pole and plane regions (Page et al. 2003a; Hinshaw et al. 2003a), and non-Gaussianity in 
the noise maps due to the 1/ / striping. Including these effects would increase our estimate of the 
power spectrum uncertainties and improve our estimate of Xeff- Our next data release will include 
the corrections and errors associated with the beam asymmetries. The features in the measured 
power spectrum could be due to underlying features in the primordial power spectrum (see §5 of 
Peiris et al. (2003)), but we do not yet attach cosmological significance to them. 

Table 1 lists the best fit parameters using the WMAP data alone for this model and Figure 
(3) shows the marginalized probabilities for each of the basic parameters in the model. The values 
in the second column of Table 1 (and the subsequent parameter tables) are expectation values for 
the marginalized distribution of each parameter and the errors are the 68% confidence interval. 
The values in the third column are the values at the peak of the likelihood function. Since we are 
projecting a high dimensional likelihood function, the peak of the likelihood is not the same as the 
expectation value of a parameter. Most of the basic parameters are remarkably well determined 
within the context of this model. Our most significant parameter degeneracy (see Figure 5) is a 
degeneracy between n s and r. The TE data favors r ~ 0.17 (Kogut et al. 2003); on the other 
hand, the low value of the quadrupole (see Figure 1 and §7) and the relatively low amplitude of 
fluctuations for I < 10 disfavors high r as reionization produces additional large scale anisotropies. 
Because of the combination of these two effects, the likelihood surface is quite flat at its peak: the 
likelihood changes by only 0.05 as r changes from 0.11 — 0.19. This particular shape depends upon 
the assumed form of the power spectrum: in §5.2, we show that models with a scale-dependent 
spectral index have a narrower likelihood function that is more centered around r = 0.17. 

Since the WMAP data allows us to accurately determine many of the basic cosmological 
parameters, we can now infer a number of important derived quantities to very high accuracy; we 
do this by computing these quantities for each model in the MCMC and use the chain to determine 
their expectation values and uncertainties. 
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Table 1. Power Law ACDM Model Parameters- WMAP Data Only 



Parameter 




Mean (68% confidence range) 


Maximum Likelihood 


Baryon Density 


n b h 2 


0.024 ±0.001 


0.023 


Matter Density 




0.14 ±0.02 


0.13 


Hubble Constant 


h 


0.72 ±0.05 


0.68 


Amplitude 


A 


0.9 ±0.1 


0.78 


Optical Depth 


T 


1 fifi +0 - 076 
u.idd_ 071 


0.10 


Spectral Index 


n s 


0.99 ±0.04 


0.97 








1431/1342 



a Fit to WMAP data only 
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Table 2. Derived Cosmological Parameters 

Parameter Mean (68% confidence range) 

Amplitude of Galaxy Fluctuations 
Characteristic Amplitude of Velocity Fluctuations 
Baryon Density /Critical Density 
Matter Density /Critical Density 
Age of the Universe 
Redshift of Reionization b 
Redshift at Decoupling 
Age of the Universe at Decoupling 
Thickness of Surface of Last Scatter 
Thickness of Surface of Last Scatter 
Redshift at Matter /Radiation Equality 
Sound Horizon at Decoupling 

Angular Diameter Distance to the Decoupling Surface 
Acoustic Angular Scale 
Current Density of Baryons 
Baryon/Photon Ratio 

a Fit to the WMAP data only 
b Assumes ionization fraction, x e = 1 
c I a = Trd c /r s 



a 8 = 0.9 ±0.1 
a 8 n^ = 0.44 ±0.10 
Q b = 0.047 ± 0.006 
Q m = 0.29 ±0.07 
t = 13.4 ± 0.3 Gyr 
z r = 17 ± 5 
z dec = W88+1 
t dec = 372 ± 14 kyr 
Az dec = 194 ± 2 
At dec = 115 ±5 kyr 

Z e q — °^ J ^-392 

r s = 144 ± 4 Mpc 
d A = 13.7 ±0.5 Gpc 
i A = 299 ± 2 

n b = (2.7 ±0.1) x 10~ 7 cm" 3 
V = (6.518;$) x io- 10 
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Table 2 lists cosmological parameters based on fitting a power law (PL) CDM model to the 
WMAP data only. The parameters i^ec an d Zd ec are determined by using the CMBFAST code 
(Seljak & Zaldarriaga 1996) to compute the redshift of the CMB "photosphere" (the peak in the 
photon visibility function). We determine the thickness of the decoupling surface by measuring 
Azdec and At^ ec , the full-width at half maximum of the visibility function. The age of the Universe 
is derived by integrating the Friedmann equation, and a$ (the linear theory predictions for the 
amplitude of fluctuations within 8 Mpc/h spheres) from the linear matter power spectrum at z = 
is computed by CMBFAST. 

4. COMPARSION WITH ASTRONOMICAL PREDICTIONS 

In this section, we compare the predictions of the best fit power law ACDM model to other 
cosmological observations. We also list in Table 10 the best fit model to the full data set: a ACDM 
model with a running spectral index (see §5.2). In particular we consider determinations of the 
local expansion rate (i.e. the Hubble constant), the amplitude of fluctuations on galaxy scales, the 
baryon abundance, ages of the oldest stars, large scale structure data and supernova la data. We 
also consider if our determination of the reionization redshift is consistent with the prediction for 
structure formation in our best fit Universe and with recent models of reionization. In §5 and 6, we 
add some of these data sets to the WMAP data to better constrain parameters and cosmological 
models. 

4.1. Hubble Constant 

CMB observations do not directly measure the local expansion rate of the Universe rather they 
measure the conformal distance to the decoupling surface and the matter-radiation ratio through 
the amplitude of the early Integrated Sachs Wolfe (ISW) contribution relative to the height of the 
first peak. For our power law ACDM model, this is enough information to "predict" the local 
expansion rate. Thus, local Hubble constant measurements are an important test of our basic 
model. 

The Hubble Key Project (Freedman et al. 2001) has carried out an extensive program of 
using Cepheids to calibrate several different secondary distance indicators (Type la supernovae, 
Tully-Fisher, Type II supernovae, and surface brightness fluctuations). With a distance modulus 
of 18.5 for the LMC, their combined estimate for the Hubble constant is Hq = 72 ± 3(stat.) ± 
7(systematic) km/s/Mpc. The agreement between the HST Key Project value and our value, 
h = 0.72 ± 0.05, is striking, given that the two methods rely on different observables, different 
underlying physics, and different model assumptions. 

As we will show in §6, models with equation of state for the dark energy very different from a 
cosmological constant (i.e., w = —1) only fit the WMAP data if the Hubble constant is much smaller 
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than the Hubble Key Project value. An independent determination of the Hubble constant that 
makes different assumptions than the traditional distance ladder can be obtained by combining 
Sunyaev-Zel'dovich and X-ray flux measurements of clusters of galaxies, under the assumption 
of sphericity for the density and temperature profile of clusters. This method is sensitive to the 
Hubble constant at intermediate redshifts (z ~ 0.5), rather than in the nearby universe. Reese et al. 
(2002), Jones et al. (2001), and Mason et al. (2001) have obtained values for the Hubble constant 
systematically smaller than, the Hubble Key Project and WMAP ACDM model determinations, 
but all consistent at the la level. Table (3) summarizes recent Hubble constant determinations and 
compares them with the WMAP ACDM model value. 

4.2. Amplitude of Fluctuations 

The overall amplitude of fluctuations on large-scale structure scales has been recently de- 
termined from weak lensing surveys, clusters number counts and peculiar velocities from galaxy 
surveys. Weak lensing surveys and peculiar velocity measurements are most sensitive to the combi- 
nation ogf^ 6 ! cluster abundance at low redshift is sensitive to a very similar parameter combination 
crgS^ 5 , but counts of high redshift clusters can break the degeneracy. 

4-2.1. Weak Lensing 

Weak lensing directly probes the amplitude of mass fluctuations along the line of sight to 
the background galaxies. Once the redshift distribution of the background galaxies is known, this 
technique directly probes gravitational potential fluctuations, and therefore can be easily compared 
with our CMB model predictions for the amplitude of dark matter fluctuations. Several groups 
have reported weak shear measurements within the past year (see Table 4 and Van Waerbeke et al. 
(2002a) for recent review): while there is significant scatter in the reported amplitude, the best fit 
model to the WMAP data lies in the middle of the reported range. As these shear measurements 
continue to improve, the combination of WMAP observations and lensing measurements will be a 
powerful probe of cosmological models. 

4-2.2. Galaxy velocity fields 

The galaxy velocity fields are another important probe of the large scale distribution of matter. 
The Willick & Strauss (1998) analysis of the Mark III velocity fields and the IRAS redshift survey 
yields (3 IRAS = 0.50 ±0.04. IRAS galaxies are less clustered than optically selected galaxies; Fisher 
et al. (1994) find al RAS = 0.69±0.04 implying <j^ ass ^ = 0.345±0.05, consistent with our ACDM 
model value of 0.44 ± 0.10. 
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Table 3. Recent Hubble Constant Determinations 



Method Mean (68% confidence range) Reference 



Hubble Key Project 


72 ± 3 ± 7 


Freedman et al. (2001) 


SZE + X-ray 


60 ± 


Reese et al. (2002) 




66tn ± 15 


Mason et al. (2001) 


WMAP PL ACDM model 


72 ±5 


§3 



Table 4. Amplitude of Fluctuations, as 



Method Mean (68% confidence range) Reference 



PL ACDM + WMAP 


0.9 ±0.1 


§3 


Weak Lensing a ' 6 


0.72 ±0.18 


Brown et al. (2002) 






Hoekstra et al. (2002) 




o.69l8:i§ 


Jarvis et al. (2002) 




0.96 ±0.12 


Bacon et al. (2002) 




0.92 ±0.2 


Refregier et al. (2002) 




0.98 ±0.12 


Van Waerbeke et al. (2002b) 


Galaxy Velocity Fields' 3 


0.73 ±0.1 


Willick & Strauss (1998) 


CBI SZ detection 


1.04±0.12 c 


Komatsu & Seljak (2002) 


High redshift clusters' 3 


0.95 ±0.1 


Bahcall k Bode (2002) 



a Since most weak lensing papers report 95% confidence limits in their papers, the table lists the 
95% confidence limit for these experiments. 

b All of the as measurements have been normalized to £l rn = 0.287, the best fit value for a fit to 
the WMAP data only. 

c 95% confidence limit 
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4-2.3. Cluster Number Counts 

Our best fit to the WMAP data is cr 8 l^ 5 = 0.48 ± 0.12. Bahcall et al. (2002b) recent study 
of the mass function of 300 clusters at redshifts 0.1 < z < 0.2 in the early SDSS data release 
yields ogl^ 5 = 0.33 ± 0.03. This difference may reflect the sensitivity of the cluster measurements 
to the conversion of cluster richness to mass. Observations of the mass function of high redshift 
clusters break the degeneracy between as and £l m . The recent Bahcall & Bode (2002) analysis of 
the abundance of massive clusters at z = 0.5 — 0.8 yields as = 0.95 ±0.1 for Vt m = 0.25. Other 
cluster analysis yield different values: Borgani et al. (2001) best fit values for a large sample of 
X-ray clusters are as = 0.66^005 and Q. m = 0.35^q'}q. On the other hand, Reiprich & Bohringer 
(2002) find very different values: a 8 = 0.96^;^ and 9, m = 0.12l°;^. Pierpaoli et al. (2002) discuss 
the wide range of values that different X-ray analyses find for as- With the larger REFLEX sample, 
Schuecker et al. (2003a) find a 8 = 0.71 ll^ and n m = 0.3411°;°^ -osnv where the second 

set of errors include the systematic uncertainties. The best fit WMAP values lie in the middle of 
the relevant range. 

Measurements of the contribution to the CMB power spectrum on small scales from the 
Sunyaev-Zel'dovich effect also probe the number density of high redshift clusters. The recent 
CBI detection of excess fluctuations (Mason et al. 2001; Bond et al. 2002) at £ > 1500 implies 
as = 1.04 ± 0.12 (Komatsu & Seljak 2002), if the signal is due to the Sunyaev-Zel'dovich effect. 



4.3. Baryon Abundance 

Both the amplitude of the acoustic peaks in the CMB spectrum (Bond & Efstathiou 1984) and 
the primordial abundance of Deuterium (Boesgaard & Steigman 1985) are sensitive functions of the 
cosmological baryon density. Since the height and position of the acoustic peaks depend upon the 
properties of the cosmic plasma 372,000 years after the Big Bang and the Deuterium abundance 
depends on physics only three minutes after the Big Bang, comparing the baryon density constraints 
inferred from these two different probes provides an important test of the Big Bang model. The best 
fit baryon abundance based on WMAP data only for the PL LCDM model, VL b h 2 = 0.024 ± 0.001, 
implies a baryon/photon ratio of r/ = (6.5^3) x 10~ 10 . For this abundance, standard big bang 
nucleosynthesis (Buries et al. 2001) implies a primordial Deuterium abundance relative to Hydrogen: 
[D]/[H] = 2.37to;2? x 10" 5 . As it will be clear from §5 and 6, the best fit VL b h 2 value for our fits is 
relatively insensitive to cosmological model and data set combination as it depends primarily on the 
ratio of the first to second peak heights (Page et al. 2003b). For the running spectral index model 
discussed in §5.2, the best fit baryon abundance, = 0.0224 ± 0.0009, implies a primordial 

[D]/[H] = 2.621^ x 10- 5 . 

How does the primordial Deuterium abundance inferred from CMB compare with that observed 
from the ISM? Galactic chemical evolution destroys Deuterium because the Deuterium nucleus is 
relatively fragile and is easily destroyed in stars. Thus, measurements of the Deuterium abundance 
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within the Galaxy are usually treated as lower limits on the primordial abundance (Epstein et al. 
1976). Local measurements of D and H absorption find [D/H] abundance near 1.5 x 10~ 5 , while 
more distant measurements by IMAP and FUSE find significant variation in Deuterium abundances 
suggesting a complex Galactic chemical history (Jenkins et al. 1999; Sonneborn et al. 2000; Moos 
et al. 2002). 

Observations of Lyman a clouds reduce the need to correct the Deuterium abundance for stellar 
processing as these systems have low (but non-zero) metal abundances. These observations require 
identifying gas systems that do not have serious interference from the Lyman a forest. The Kirkman 
et al. (2003) analysis of QSO HS 243+3057 yields a D/H ratio of 2.42±g;|| x 10" 5 . They combine this 
measurement with four other D/H measurements (Q0130-4021: D/H< 6.8 x 10~ 5 , Q1009+2956: 
3.98±0.70x 10~ 5 , PKS 1937-1009: 3.25±0.28 x 10~ 5 , and QSO HS0105+1619: 2.5±0.25 x 10~ 5 ), to 
obtain their current best D/H ratio: 2.78±q;|| x 10~ 5 implying Q b h 2 = 0.0214 ± 0.0020. D'Odorico 
et al. (2001) find 2.24 ± 0.67 x 10" 5 from their observations of Q0347-3819 (although a reanalysis 
of the system by Levshakov et al. (2003) finds a higher D/H value: 3.75 ± 0.25. Pettini & Bowen 
(2001) report a D/H abundance of 1.65 ±0.35 x 10~ 5 from STIS measurements of QSO 2206-199, a 
low metallicity (Z ~ 1 /200) Damped Lyman a system. The WMAP value lies between the Pettini 
& Bowen (2001) estimate from DLAs, n b h 2 = 0.025±0.001, and the Kirkman et al. (2003) estimate 
of r^/i 2 = 0.0214 ± 0.0020 The remarkable agreement between the baryon density inferred from 
D/H values and our measurements is an important triumph for the basic Big Bang model. 



4.4. Cosmic Ages 

The age of the Universe based on the best fit to WMAP data only, to = 13.4±0.3 Gyr. However, 
the addition of other data sets (see §5) implies a lower matter density and a slightly larger age. 
The best fit age for the power law model based on a combination of WMAP , 2dFGRSand Lyman 
a forest data is to = 13.6 ± 0.2Gyr. The best fit age for the same data set for the running index 
model of §5.2 is to = 13.7 ± 0.2Gyr. (See Hu et al. 2001 Hu (Fukugita) and Knox, Christensen &z 
Skordis (2001) for discussions of using CMB data to determine cosmological ages.) 

A lower limit to the age of the universe can independently be obtained from dating the oldest 
stellar populations. This has been done traditionally by dating the oldest stars in the Milky Way 
(see e.g., Chaboyer (1998); Jimenez (1999)). For this program, globular clusters are an excellent 
laboratory for constraining the age of the universe: each cluster has a chemically homogeneous 
population of stars all born nearly simultaneously. The main uncertainty in the age determination 
comes from the poorly known distance (Chaboyer 1995). Well-understood stellar populations are 
useful tools for constraining cluster distances: Renzini et al. (1996) used the white dwarf sequence to 
obtain an age of 14.5 ± 1.5 Gyr for NGC 6752. Jimenez et al. (1996). using a distance-independent 
method determined the age of the oldest globular clusters to be 13.5 ± 2 Gyr. Using the luminosity 
function method, Jimenez & Padoan (1998) found an age of 12.5 ± 1.0 Gyr for M55. This method 
gives a joint constraint on the distance and the age of the globular cluster. Other groups find 
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Table 5. Measured ratio of Deuterium to Hydrogen 



Quasar [D]/[H] Reference 



Q0130-403 


< 6.8 x 10~ 5 


Kirkman et al. (2000) 


PKS 1937-1009 


3.25 ±0.3 x 10~ 5 


Buries & Tytler (1998a) 


Q1009+299 


4.0 ± 0.65 x 10~ 5 


Buries & Tytler (1998b) 


HS0105+1619 


2.5 ±0.25 x 10~ 5 


O'Meara et al. (2001) 


Q2206-199 


1.65 ±0.35 x 10~ 5 


Pettini k Bowen (2001) 


Q0347-383 


3.75 ±0.25 x 10~ 5 


Levshakov et al. (2003) 


Q1234+3047 


2.42±g;|| x 10~ 5 


Kirkman et al. (2003) 



Table 6. Cosmic Age 



Method 


Age 


WMAP data (ACDM) 


13.4 ± 0.3 Gyr 


WMAPext± LSS 


13.7±0.2Gyr 


Globular Cluster Ages 


> 11 - 16 Gyr 


White Dwarf 


> 12.7 ±0.7 Gyr 


OGLEGC-17 


> 10.4 - 12.8 Gyr 


Radioactive dating 


> 9.5 - 20 Gyr 
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consistent ages: Gratton et al. (1997) estimate an age of 11.8jl 2 ' 5 Gyr for the oldest Galactic 
globulars; VandenBerg et al. (2002) estimates an age of ~ 13.5 Gyr for M92. Chaboyer & Krauss 
(2003) review the globular cluster analysis and quote a best fit age of 13.4 Gyr. 

Observations of eclipsing double line spectroscopic binaries enable globular cluster age deter- 
minations that avoid the considerable uncertainty associated with the globular cluster distance 
scale (Paczynski 1997). Thompson et al. (2001) were able to obtain a high precision mass estimate 
for the detached double line spectroscopic binary, OGLEGC-17 in uo— Cen. Using the age/turnoff 
mass relationship, the Kaluzny et al. (2002) analysis of this system yielded an age for this binary 
of 11.8 ± 0.6 Gyr. Chaboyer & Krauss (2002) re-analysis of the age/turnoff mass relationship for 
this system yields a similar age estimate: 11.1 ± 0.67 Gyr. The WMAP determination of the age 
of the universe implies that globular clusters form within 2 Gyr after the Big Bang, a reasonable 
estimate that is consistent with structure formation in the ACDM cosmology. White dwarf dating 
provides an alternative approach to the traditional studies of the main sequence turn-off. Richer 
et al. (2002) and Hansen et al. (2002) find an age for the globular cluster M4 of 12.7 ± 0.7 Gyrs (2 
a errors, ±0.35 at the 1 a level assuming Gaussian errors) using the white dwarfs cooling sequence 
method. These results, which yield an age close to the cosmological age, are potentially very useful: 
further tests of the assumptions of the white dwarf age dating method will clarify its systematic 
uncertainties. 

Observations of nearby halo stars enable astronomers to obtain spectra of various radio- 
isotopes. By measuring isotopic ratios, they infer stellar ages that are independent of much of 
the physics that determines main sequence turn-off (see Thielemann et al. (2002) for a recent re- 
view). These studies yield stellar ages consistent with both the globular cluster ages and the ages 
in our best fit models. Clayton (1988) using a range of chemical evolution models for the Galaxy 
finds ages between 12 - 20 Gyr. Schatz et al. (2002) study Thorium and Uranium in CS 31082-001 
and estimate an age of 15.5 ± 3.2 Gyr for the r-process elements in the star. Other groups find 
similar estimates: the Cayrel et al. (2001) analysis of U-238 in the old halo star CS 31082-001 yields 
an age of 12.5 ± 3 Gyr, while Hill et al. (2002) find an age of 14.0 ± 2.4 Gyr. Studies of other old 
halo stars yield similar estimates: Cowan et al. (1999) two stars CS 22892-052 and HD115444 find 
15.6 ±4.6 Gyr. 

Table 6 summarizes the lower limits on the age of the universe from various astronomical 
measurements. While the errors on these measurements remain too large to effectively constrain 
parameters, they provide an important consistency check on our basic cosmological model. 

4.5. Large Scale Structure 

The large scale structure observations and the Lyman a forest data complement the CMB 
measurements by measuring similar physical scales at very different epochs. The WMAP angular 
power spectrum has the smallest uncertainties near I ~ 300, which correspond to wavenumbers 
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k ~ 0.02 Mpc 1 . With the AC BAR results, our CMB data set extends to I ~ 1800, corresponding 
to k ~ 0.1 Mpc" 1 . If we assume that gravity is the primary force determining the large-scale 
distribution of matter and that galaxies trace mass at least on large scales, then we can directly 
compare our best fit A CDM model (with parameters fit to the WMAP data) to observations of large 
scale distribution of galaxies. There are currently two major ongoing large scale structure surveys: 
the Anglo-Australian Telescope two degree field Galaxy Redshift Survey (2dFGRS) (Colless et al. 
2001), and the Sloan Digital Sky Survey 14 (SDSS). Large scale structure data sets are a powerful 
tool for breaking many of the parameter degeneracies associated with CMB data. In §5, we make 
extensive use of the 2dFGRS data set. 

Figure 6 shows that the ACDM model obtained from the WMAP data alone is an acceptable 
fit to the 2dFGRS power spectrum. The best fit has j3 = 0.45 consistent with Peacock et al. (2001) 
measured value of (3 = 0.43 ± 0.07. 

The Lyman a forest observations are an important complement to CMB observations since they 
probe the linear matter power spectrum at z = 2— 3 (Croft et al. 1998, 2002). These observations are 
sensitive to small length scales, inaccessible to CMB experiments. Unfortunately, the relationship 
between the measured flux power spectrum and the linear power spectrum is complex (Gnedin & 
Hamilton 2002; Croft et al. 2002) and needs to be calibrated by numerical simulations. In Verde 
et al. (2003), we describe our methodology for incorporating the Lyman a forest data into our 
likelihood approach. Figure 6 compares the predicted power spectra for the best fit ACDM model 
to the linear power spectra inferred by Gnedin & Hamilton (2002) and by Croft et al. (2002). 

4.6. Supernova Data 

Over the past decade, Type la supernovae have emerged as important cosmological probes. 
Once supernova light curves have been corrected using the correlation between decline rate and 
luminosity (Phillips 1993; Riess et al. 1995) they appear to be remarkably good standard candles. 
Systematic studies by the supernova cosmology project (Perlmutter et al. 1999) and by the high 
z supernova search team (Riess et al. 1998) provide evidence for an accelerating universe. The 
combination of the large scale structure, CMB and supernova data provide strong evidence for a 
flat universe dominated by a cosmological constant (Bahcall et al. 1999). Since the supernova data 
probes the luminosity distance versus redshift relationship at moderate redshift z < 2 and the CMB 
data probes the angular diameter distance relationship to high redshift (z ~ 1089) , the two data 
sets are complementary. The supernova constraint on cosmological parameters are consistent with 
the ACDM WMAP model. As we will see in the discussion of non-flat models and quintessence 
models, the SNIa likelihood surface in the Q, m — Q\ and in the Q m — w planes provides useful 
additional constraints on cosmological parameters. 



www.sdss.org 
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4.7. 



Reionization & Small Scale Power 



The WMAP detection of reionization (Kogut et al. 2003) implies the existence of an early 
generation of stars able to reionize the Universe at z ~ 20. Is this early star formation compatible 
with our best fit ACDM cosmological model? We can evaluate this effect by first computing the 
fraction of collapsed objects, /dm, at a given redshift: 



where <3?(M, z) is the Sheth & Tormen (1999) mass function. The first stars correspond to extremely 
rare fluctuations of the overdensity field: Eq. (2) is very sensitive to the tail of the mass function. 
Thus the very small change in the minimum mass needed for star formation results in a significant 
change in the fraction of collapsed objects. The minimum halo mass for star formation, M m i n , is 
controversial and depends on whether molecular hydrogen (H2) is available as a coolant. If the gas 
temperature is fixed to the CMB temperature, then the Jean Mass, M- 7 = 10 6 M Q . If molecular 
hydrogen is available, then the Jeans mass before reionization is M J ' ~ 2.2 x 10^[tUb/h(uj m )] l ' 5 (l + 
z)/10 for z < 150 (Venkatesan et al. 2001). At z > 150, the electrons are thermally coupled to the 
CMB photons. However, as Haiman et al. (1997) point out, a small UV background generated by 
the first sources will dissociate H2, thus making the minimum mass much larger than the Jeans 
mass. They suggest using a minimum mass that is much higher: M^^ L (z) = 10 8 (1 + z)/10)~ 3 / 2 . 
On the other hand if the first stars generated a significant flux of X-rays (Oh 2001) then this would 
have promoted molecular hydrogen formation (Haiman et al. 2000; Venkatesan et al. 2001; Cen 
2002). Thus lowering the minimum mass back to MK 

Following Tegmark & Silk (1995) we estimate the rate of reionization by multiplying the col- 
lapse factor by an efficiency factor. A fraction of baryons in the universe, fb, falls into the non-linear 
structures. We assume fb = /dm (i-e., constant baryon/dark matter ratio). A certain fraction of 
these baryons form stars or quasars, fburm which emit UV radiation with some efficiency, fuv- Some 
of this radiation escapes into the intergalactic medium photoionizing it; however, the net number 
of ionizations per UV photons, /j on , is expected to be less than unity (due to cooling and recombi- 
nations). Finally the intergalactic medium might be clumpy, making the photoionization process 
less efficient. This effect is counted for by the clumping factor C c i ump . Thus in this approximation 
the ionization fraction is given by: x e = 3.8 x 10 5 fnetfb where f net = fburnfuv f esc f ion/ 'C dump . The 
factor 3.8 x 10 5 arises because 7.3 x 10 3 of the rest mass is released in the burning of hydrogen to 
helium and we assume the primordial helium mass fraction to be 24%. We assume fbum ^5 25%, 



Figure 7 shows the fraction of collapsed objects and the maximum ionization fraction as a 
function of redshift for our best fit WMAP ACDM model. The solid lines correspond to M m i n = 
M^^ L (z) while the dashed lines correspond to M m i n = M- 7 . The WMAP detection of reionization 
at high redshift suggests that H2 cooling likely played an important role in early star formation. 




(2) 



fesc < 50%, f uv < 50%, f lon < 90%, and 



1 < C dump < 100, thus f net < 5.6 x 10- 3 . 



Because early reionization requires the existence of small scale fluctuations, the WMAP TE 
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detection has important implications for our understanding of the nature of the dark matter. 
Barkana et al. (2001) note that the detection of reionization at z > 10 rules out warm dark matter 
as a viable candidate for the missing mass as structure forms very late in these models. Warm dark 
matter can not cluster on scales smaller than the dark matter Jeans' mass. Thus, this limit applies 
regardless of whether the minimum mass is M HRL or MK 



5. COMBINING DATA SETS 

In this section, we combine the WMAP data with other CMB experiments that probe smaller 
angular scales (ACBAR and CBI) 15 and with astronomical measurements of the power spectrum 
(the 2dFGRS and Lyman a forest). We begin by exploring how including these data sets affects our 
best fit power law ACDM model parameters (§5.1). The addition of data sets that probe smaller 
scales systematically pulls down the amplitude of the fluctuations in the best fit model. This 
motivates our exploration of an extension of the power law model, a model where the primordial 
power spectrum of scalar density fluctuations is fit by a running spectral index (Kosowsky & Turner 
1995): 

/ i \ n s (feo)+(l/2)dn s /<ilnfeln(fe/A;o) 

P(k)=P(k ) -) , (3) 



where we fix the scalar spectral index and slope at k® = 0.05Mpc _1 . Note that this definition of the 
running index matches the definition used in Hannestad et al. (2002) analysis of running spectral 
index models and differs by a factor of 2 from the Kosowsky & Turner (1995) definition. As in the 
scale independent case, we define 

We explicitly assume that d 2 n s /dlnk 2 = 0, so that 

„.(*)_„.<„,)+ *L ,„(£). (5) 

In §5.2, we show that the running spectral index model is a better fit than the pure power 
law model to the combination of WMAP and other data sets. Peiris et al. (2003) explores the 
implications of this running spectral index for inflation. 



5.1. Power Law CDM Model 



The power law ACDM model is an acceptable fit to the WMAP data. While it overpredicts 
the amplitude of fluctuations on large angular scales (see §6), this deviation may be due to cosmic 



In the following sections, we refer to the combined WMAP , ACBAR and CBI data sets as WM APext. 
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variance at these large scales. Intriguingly, it also overpredicts the amplitude of fluctuations on 
small angular scales. 

Table (7) shows the best fit parameters for the power law ACDM model for different combi- 
nation of data sets. As we add more and more data on smaller scales, the best fit value for the 
amplitude of fluctuations at k = 0.05 Mpc -1 gradually drops: When we fit to the WMAP data 
alone, the best fit is 0.9 ± 0.1. When we add the CBI, ACBAR and 2dFGRS data, the best fit 
value drops to 0.8 ± 0.1. Adding the Lyman a data further reduces A to 0.75+QQ7. The best fit 
spectral index shows a similar trend: the addition of more and more small scale data drives the 
best fit spectral index to also change by nearly 1 a from its best fit value for WMAP data only: 
0.99 ± 0.04 (WMAP only) to 0.96 ± 0.02 ( WMAPext+2dFGRS+Lya). When the addition of new 
data continuously pulls a model away from its best fit value, this is often the signature of the model 
requiring a new parameter. 

5.2. Running Spectral Index ACDM Model 

Inflationary models predict that the spectral index of fluctuations should be a slowly varying 
function of scale. Peiris et al. (2003) discusses the inflationary predictions and shows that a plausible 
set of models predicts a detectable varying spectral index. There are classes of inflationary models 
that predict minimal tensor modes. This section explores this class of models. In §6.4, we explore 
a more general model that has both a running spectral index and tensor modes. 

Table 8 shows the best fit parameters for the running (RUN) spectral index model as a function 
of data set. Note that the best fit parameters for these models barely change as we add new data 
sets; however, the error bars shrink. When we include all data sets, the best fit value of the running 
of the spectral index is — O.OSIIq'q}^: fewer than 5% of the models have dn s /dlnk > 0. 

Figure 9 shows the the power spectrum as a function of scale. The figure shows the results of 
our Markov chain analysis of the combination of WMAP , CBI, ACBAR, 2dFGRS and Lyman a 
data. At each wavenumber, we compute the range of values for the power law index for all of the 
points in the Markov chain. The 68% and 95% contours at each k value are shown in Figure 9 for 
the fit to the WM APext + 2dFGRS + Lyman a data sets. 

Over the coming year, new data will significantly improve our ability to measure (or constrain) 
this running spectral index. When we complete our analysis of the EE power spectrum, the WMAP 
data will place stronger constraints on r. Because of the n s — r degeneracy, this implies a strong 
constraint on n s on large scales. The SDSS collaboration will soon release its galaxy spectrum 
and its measurements of the Lyman a forest. These observations will significantly improve our 
measurements of n s on small scales. Peiris et al. (2003) shows that the detection of a running 
spectral index and particularly the detection of a spectral index that varies from n s > 1 on large 
scales to n s < 1 on small scales would severely constrain inflationary models. 



Table 7. Best Fit Parameters: Power Law A CDM 



WMAP WMAPext 16 a WMAPext+2dFGRS WMAPext+ 2dFGRS+ Lyman a 



A 


0.9 ±0.1 


0.8 ±0.1 


0.8 ±0.1 


n 7 c:+0.08 
u - '°-0.07 


n s 


0.99 ± 0.04 


0.97 ±0.03 


0.97 ±0.03 


0.96 ±0.02 


T 


0.166lg;g» 


n 143+0-071 

U.140_q 062 


o.i48i8:8™ 


n 1 1 7 +0.057 
U- 11 ' -0.053 


h 


0.72 ± 0.05 


0.73 ±0.05 


0.73 ± 0.03 


0.72 ± 0.03 


Quill 


0.14 ±0.02 


0.13 ±0.01 


0.134 ±0.006 


0.133 ±0.006 


n b h 2 


0.024 ± 0.001 


0.023 ± 0.001 


0.023 ±0.001 


0.0226 ± 0.0008 


Xlff/ U 


1429/1341 


1440/1352 


1468/1381 


b 



a WMAP ±CBI±ACBAR 

b Since the Lyman a data points are correlated, we do not quote an effective \ 2 f° r the combined 
likelihood including Lyman a data (see Verde et al. (2003)). 



Table 8. Best Fit Parameters for the Running Spectral Index ACDM Model 



WMAP WMAPext WMAPext+2dFGRS WMAPcxt+ 2dFGRS+ Lyman a 



A 


0.92 + 0.12 


0.9 ±0.1 


0.84 ±0.09 


n co+o.09 


n 3 


0.93l»; ' 7 7 


0.91 ±0.06 


n qo+o.04 


0.93 ±0.03 


dris/dln k 


-0.047 ± 0.04 


-0.055 ±0.038 


-0 031 + 023 
u.uoi_ 025 


-0 031 + 016 


T 


0.20 ±0.07 


0.20 ±0.07 


0.17 ±0.06 


0.17 ±0.06 


h 


0.70 ± 0.05 


0.71 ± 0.06 


0.71 ±0.04 


71 + ' 04 
u - ' J^-o.oa 




0.14 + 0.02 


0.14 ± 0.01 


0.136 ±0.009 


n i qc+o.008 
u.iod_ 009 


Q b h 2 


0.023 ±0.002 


0.022 ± 0.001 


0.022 ±0.001 


0.0224 + 0.0009 




1431/1342 


1437/1350 


1465/1380 





a Since the Lyman a data points are correlated, we do not quote Xeff f° r the combined likelihood including 
Lyman a data (see Verde et al. (2003)). 
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The running spectral index model predicts a significantly lower amplitude of fluctuations on 
small scales than the standard ACDM model (see figure 9). This suppression of small scale power 
has several important astronomical implications: (a) the reduction in small scale power makes it 
more difficult to reionize the universe unless Hi cooling enables mass dark halos to collapse and form 
galaxies (see §4.7 and Figure 10); (b) a reduction in the small scale power reduces the amount of 
substructure within galactic halos (Zentner & Bullock 2002) (c) since small objects form later, their 
dark matter halos will be less concentrated as there is a monotonic relationship between collapse 
time and halo central concentration (Navarro et al. 1997; Eke et al. 2001; Zentner & Bullock 2002; 
Wechsler et al. 2002; Huffenberger & Seljak 2003). The reduction in the amount of substructure will 
also reduce angular momentum transport between dark matter and baryons and will also reduce 
the rate of disk destruction through infall (Toth & Ostriker 1992). We suspect that our proposed 
modification of the primordial power spectrum will resolve many of the long-standing problems of 
the CDM model on small scales (see Moore (1994) and Spergel h Steinhardt (2000) for discussions 
of the failings of the power law A CDM model on galaxy scales). 

6. BEYOND THE ACDM MODEL 

In this section, we consider various extensions to the ACDM model. In §6.1, we consider dark 
energy models with a constant equation of state. In §6.2, we consider non-flat models. In §6.3, we 
consider models with a massive light neutrino. In §6.4, we include tensor modes. 

In this section of the paper, we combine the WMAP data with external data sets so that we 
can break degeneracies and obtain significant constraints on the various extensions of our standard 
cosmological model. 

6.1. Dark Energy 

The properties of the dark energy, the dominant component in our universe today, is a mystery. 
The most popular alternative to the cosmological constant is quintessence. Wetterich (1988), Ratra 
& Peebles (1988) and Peebles &; Ratra (1988) suggest that a rolling scalar field could produce a 
time-variable dark energy term, which leave a characteristic imprint on the CMB and on large 
scale structure (Caldwell et al. 1998). In these quintessence models, the dark energy properties are 
quantified by the equation of state of the dark energy: w = p/p, where p and p are the pressure 
and the density of the dark energy. A cosmological constant has an equation of state, w = — 1. 

Since the space of possible models is quite large, we only consider models with a constant 
equation of state. We now increase our model space so that we have 7 parameters in the cosmological 
model (A, n s , h, J7 m , f^, r and w). We analyze the data using two approaches: (a) we begin by 
restricting our analysis to w > —1 motivated by the difficulties in constructing stable models with 
w < —I (Carroll et al. 2003) and (b) relax this constraint and consider models that violate the 
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weak energy condition (Schuecker et al. 2003b). Further analysis is needed for models where w 
and the quintessence sound speed are a function of time (Dedeo et al. 2003). The addition of a new 
parameter introduces a new degeneracy between £l m , h, and w that can not be broken by CMB 
data alone (Huey et al. 1999; Verde et al. 2003): models with the same values of rt m h 2 , Vt^h? and 
first peak position have nearly identical angular power spectra. 

For example, a model with £l m = 0.47, w = —1/2 and h = 0.57 has a nearly identical angular 
power spectrum to our ACDM model. Note, however, that this Hubble Constant value differs by 
2(7 from the HST Key Project value and the predicted shape of the power spectrum is a poor fit 
to the 2dFGRS observations. This model is also a worse fit to the supernova angular diameter 
distance relation. 

We consider four different combinations of astronomical data sets: (a) WMAPext data com- 
bined with the supernova observations; (b) WMAPext data combined with HST data; (c) WMAPext 
data combined with the 2dFGRS large scale structure data; (d) all data sets combined. 

The CMB peak positions constrain the conformal distance to the decoupling surface. The 
amplitude of the early ISW signal determines the matter density, Q m h 2 . The combination of these 
two measurements strongly constrains £l(w) and h(w) (see Figures 11 and 12). The HST Key 
Project measurement of H$ agrees with the inferred CMB value if w = — 1. As w increases, the 
best fit Hq value for the CMB drops below the Key Project value. Our joint analysis of CMB + 
HST Key Project data implies that w < —0.5 (95% confidence interval). If future observations 
can reduce the uncertainties associated with the distance to the LMC, the Hq measurements could 
place significantly stronger limits on w. Figures 11 and 12 show that the combination of either 
CMB+supernova data or CMB+large scale structure data place similar limits on dark energy 
properties. For our combined data set, we marginalize over all other parameters and find that 
w < —0.78 (95% CL) when we impose the prior that w > —1. If we drop this prior, then all of the 
combined data sets appear to favor a model where the properties of the dark energy are close to 
the predicted properties of a cosmological constant (w = —0.98 ± 0.12). 



6.2. Non-Flat Models 

The position of the first peak constrains the universe to be nearly flat (Kamionkowski et al. 

— 1/2 

1994); low density models with = have their first peak position at I ~ 200Q m . However, if 
we allow for the possibility that the universe is non-flat and there is a cosmological constant, then 
there is a geometric degeneracy (Efstathiou & Bond 1999): along a line in £l m — Qa space, there 
is a set of models with nearly identical angular power spectra. While the allowed range of Q, tot is 
relatively small, there is a wide range in VL m values compatible with the CMB data in a non-flat 
universe. 

If we place no priors on cosmological parameters, then there is a model with £l\ = consistent 
with the WMAP data (Ax 2 = 6.6 relative to the flat model). However, the cosmological parameters 
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for this model {Hq = 32.5 km/s/Mpc, and £l t ot = 1-28) are violently inconsistent with a host of 
astronomical measurements. The flat £l m = 1, A = standard CDM model is inconsistent with the 
WMAP data at more than the 5a level. 

If we include a weak prior on the Hubble Constant, Hq > 50km/s/Mpc, then this is sufficient 
to constrain 0.98 < Qtot < 1-08 (95% confidence interval). Combining the WMAPextdata with 
supernova measurements of the angular diameter distance relationship (see figure 13) we obtain 
0.98 < £ltot < 1-06. This confidence interval does not require a prior on h. If we further include 
the HST Key Project measurement of Hq as a prior, then the limits on Qq improve slightly: 
Qtot < 1.02 ± 0.02 Figure 13 shows the two dimensional likelihood surface for various combinations 
of the data. 



6.3. Massive Neutrinos 

Copious numbers of neutrinos were produced in the early universe. If these neutrinos have 
non-negligible mass they can make a non-trivial contribution to the total energy density of the 
universe during both matter and radiation domination. During matter domination, the massive 
neutrinos cluster on very large scales but free-stream out of smaller scale fluctuations. This free- 
streaming changes the shape of the matter power spectrum (Hu et al. 1998) and most importantly, 
suppresses the amplitude of fluctuations. Since we can normalize the amplitude of fluctuations to 
the WMAP data, the amplitude of fluctuations in the 2dFGRS data places significant limits on 
neutrino properties. 

The contribution of neutrinos to the energy density of the universe depends upon the sum of 
the mass of the light neutrino species: 

Oh 2 - ^ iTni (6) 
vh ~94.0eV [ ) 

Note that the sum only includes neutrino species light enough to decouple while still relativistic. 

Experiments that probe neutrino propagation from source to detector are sensitive not to the 
neutrino mass but to the square mass difference between different neutrino mass eiginstates. Solar 
neutrino experiments (Bahcall et al. 2003a) imply that the square mass difference between the 
electron and muon neutrino is ~ 10~ 9 eV. The deficit of muon neutrinos in atmospheric showers 
imply that the mass difference between muon and tau neutrinos is 10 _5 eV 2 (Kearns 2002). If 
the electron neutrino is much lighter than the tau neutrino, then the combination of these results 
imply that vtl Vt < 0.1 eV: still below the detection limits for our data-set. On the other hand, if 
m Ue ~ m„ T , then the three neutrino species can leave an observable imprint on the CMB angular 
power spectrum and the galaxy large scale structure power spectrum. In our analysis, we consider 
this latter case and assume that there are three degenerate stable light neutrino species. 

Figure 14 shows the cumulative likelihood of the combination of WMAP, CBI, ACBAR, and 
2dFGRS data as a function of the energy density in neutrinos. Based on this analysis, we conclude 
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that Q u h 2 < 0.0067 (95% confidence limit). If we add the Lyman a data, then the limit slightly 
weakens to VL u h 2 < 0.0076. For three degenerate neutrino species, this implies that m v < 0.23 eV. 
This limit is roughly a factor of two improvement over previous analyses (e.g., Elgar0y et al. (2002)) 
that had to assume strong priors on Q m and -f^o- 



6.4. Tensors 

Many models of inflation predict a significant gravity wave background. These tensor fluctua- 
tions were generated during inflation. Tensor fluctuations have their largest effects on large angular 
scales where they add in quadrature to the fluctuations generated by scalar modes. 

Here, we place limits on the amplitude of tensor modes. We define the tensor amplitude using 
the same convention as Leach et al. (2002): 

^ Ptensorik*) 

Pscalar{k*) 

where Ptensor and P sca iar are the primordial amplitude of tensor and scalar fluctuations and = 
0.002 Mpc -1 . Since we see no evidence for tensor modes in our fit, we simplify the analysis by 
assuming that the tensor spectral index satisfies the single field inflationary consistency condition: 

n t = -r/8. (8) 

This constraint reduces the number of parameters in this model to 8: A, ^h 2 , £l m h 2 , h, n s , 
dn s /dlnk, r and r. We ignore the running of n t . The addition of this new parameter does not 
improve the fit as figure (15) shows the combination of WMAPext+ 2dFGRS + Lyman a is able 
to place a limit on the tensor amplitude: r < 0.90 (95% confidence limit). As table (9) shows, this 
limit is much more stringent if we restrict the parameter space to models with either n s < 1 or 
\dn/d\nk\ = 0. 

Peiris et al. (2003) discuss the implications of our limits on tensor amplitude for inflationary 
scenarios. Using the results of this analysis, Peiris et al. (2003) shows that the inferred joint 
likelihood of n s , dn s /d\nk and r places significant constraints on inflationary models. 



Table 9. 95% Confidence Limits on Tensor/Scalar Ratio 



prior 


WMAP 


WMAPext+2dFGRS 


WMAPext+ 2dFGRS+ Lyman a 


no prior 


1.28 


1.14 


0.90 


dn s /d In k = 


0.81 


0.53 


0.43 


n s < 1 


0.47 


0.37 


0.29 
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7. INTRIGUING DISCREPANCIES 

While the ACDM model's success in fitting CMB data and a host of other astronomical data is 
truly remarkable, there remain a pair of intriguing discrepancies: on both the largest and smallest 
scales. While adding a running spectral index may resolve problems on small scales, there remains 
a possible discrepancy between predictions and observations on the largest angular scales. 

Figure 16 shows the measured angular power spectrum and the predictions of our best fit 
A— CDM model, where the data were fit to both CMB and large-scale structure data. The figure 
also shows the measured angular correlation function; the lack of any correlated signal on angular 
scales greater than 60 degrees is noteworthy. We quantify this lack of power on large scales by 
measuring a four point statistic: 

S = J [C '(9)} 2 d cos 6. (9) 

The upper cutoff and the form of this statistic were both determined a posteori in response to 
the shape of the correlation function. We evaluate the statistical significance of these discrepancies 
by doing Monte-Carlo realizations of the first 100,000 models in the Markov chains. This allows us 
to average not only over cosmic variance but also over our uncertainties in cosmological parameters. 
For our ACDM Markov chains (fit to the WMAPext+ 2dFGRS data sets), we find that only 0.7% 
of the models have lower values for the quadrupole and only 0.15% of the simulations have lower 
values of S. For the running model, we find that only 0.9% of the models have lower values for 
the quadrupole and only 0.3% of the simulations have lower values of S. The shape of the angular 
correlation function is certainly unusual for realizations of this model. 

Is this discrepancy meaningful? The low quadrupole was already clearly seen in COBE and 
was usually dismissed as due to cosmic variance (Bond et al. 1998) or foreground contamination. 
While the WMAP data reinforces the case for its low value, cosmic variance is significant on these 
large angular scales and any Gaussian field will always have unusual features. On the other hand, 
this discrepancy could be the signature of interesting new physics. 

The discovery of an accelerating universe implies that at these large scales, there is new and 
not understood physics. This new physics is usually interpreted to be dark energy or a cosmological 
constant. In either case, we would expect that the decay of fluctuations at late times produces a 
significant ISW signal. Boughn et al. (1998) argue that in a ACDM model with Q m = 0.25, there 
should be a detectable correlation between the CMB signal and tracers of large-scale structure; 
yet they were not able to detect a signal. There are alternative explanations of the accelerating 
universe, such as extra dimensional gravity theories (Deffayet et al. 2002) that do not require a 
cosmological constant and should make radically different predictions for the CMB on these angular 
scales. These predictions have not yet been calculated. 

What could generate this unusual shaped angular correlation function? As an example, we 
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compute the angular correlation function in a toy model, where the power spectrum has the form: 



where tq is the conformal distance to the surface of last scatter. This toy model simulates both 
the effects of a discrete power spectrum due to a finite universe and the effects of ringing in the 
power spectrum due to a feature in the inflaton potential (see Peiris et al. (2003) for a discussion 
of inflationary models). Figure 16 shows the angular correlation function and figure 17 show the 
TE power spectrum of the model. Note that the TE power spectrum is particularly sensitive 
to features in the matter power spectrum. Intriguingly, this toy model is a better match to the 
observed correlation function than the ACDM model and predicts a distinctive signature in the TE 
spectrum. Cornish et al. (1998) show that if the universe was finite and smaller than the volume 
within the decoupling surface, then there should be a very distinctive signal: matched circles. The 
surface of last scatter is a sphere centered around WMAP. If the universe is finite then this sphere 
must intersect itself, this leads to pairs of matched circles. These match circles provide not only 
the definitive signature of a finite universe but also should enable cosmologists to determine the 
topology of the universeCornish et al. (1998b); Weeks (1998). Should we be able to detect circles 
if the power spectrum cutoff is due to the size of the largest mode being ~ 1/ro? While there is no 
rigorous theorem relating the size of the largest mode to the diameter of the fundamental domain, 
D, analysis of both negatively curved (Cornish & Spergel 2000) and positively curved (Lehoucq 
et al. 2002) topologies suggest that D ~ (0.6 — 1)A. Thus, if the "peak" in the power spectrum 
at / = 5 corresponds to the largest mode in the domain, we should be able to detect a pattern of 
circles in the sky. 

Due to the finite size of the patch of the universe visible to WMAP (or any future satellite) , our 
ability to determine the origin and significance of this discrepancy will be limited by cosmic variance. 
However, future observations can offer some new insight into its origin. By combining the WMAP 
data with tracers of large scale structure (Boughn et al. 1998; Peiris & Spergel 2000), astronomers 
may be able to directly detect the component of the CMB fluctuations due to the ISW effect. 
WMAP's ongoing observations of large-scale microwave background polarization fluctuations will 
enable additional measurements of fluctuations at large angular scales. Since the TE observations 
are probing different regions of the sky from the TT observations, they may enlighten us on whether 
the lack of correlations on large angular scales is a statistical fluke or the signature of new physics. 



Cosmology now has a standard model: a flat universe composed of matter, baryons and vacuum 
energy with a nearly scale-invariant spectrum of primordial fluctuations. In this cosmological 
model, the properties of the universe are characterized by the density of baryons, matter and the 
expansion rate: f^, Q m , and h. For the analysis of CMB results, all of the effects of star formation 




(10) 



n=l 



8. CONCLUSIONS 
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Table 10. Basic and Derived Cosmological Parameters: Running Spectral Index Model 3, 



Mean and 68% Confidence Errors 



Amplitude of fluctuations 


A = 0.8318:" 


Spectral Index at k = 0.05 Mpc -1 


n s = 0.93 ± 0.03 


Derivative of Spectral Index 


dn s /d\nk = -0.03ll{j;{£g 


Hubble Constant 


u — n vi+0.04 
a — u.i i_ 03 


Baryon Density 


tt b h 2 = 0.0224 ± 0.0009 


Matter Density 


n U2 _ n -|or;+0.008 


Optical Depth 


r = 0.17 ±0.06 


Matter Power Spectrum Normalization 


cr 8 = 0.84 ± 0.04 


Characteristic Amplitude of Velocity Fluctuations 


^ 0O.6 n qo+0.04 

a 8 n^ = 0.38_ 05 


Baryon Density/Critical Density 


Q b = 0.044 ± 0.004 


Matter Density/Critical Density 


Q rn = 0.27 ±0.04 


Age of the Universe 


t = 13.7 ±0.2 Gyr 


Reionization Redshift b 


z r = 17 ± 4 


Decoupling Redshift 


z dec = 1089 ± 1 


Age of the Universe at Decoupling 


t dec = 3791? kyr 


Thickness of Surface of Last Scatter 


Az dec = 195 ± 2 


Thickness of Surface of Last Scatter 


At dec = 118l| kyr 


Redshift of Matter/Radiation Equality 


%eq = 3233 210 


Sound Horizon at Decoupling 


r s = 147 ± 2 Mpc 


Angular Size Distance to the Decoupling Surface 


d A = 14.0±g;| Gpc 


Acoustic Angular Scale 


t A = 301 ± 1 


Current Density of Baryons 


n b = (2.5 ±0.1) x lO" 7 cm" 3 


Baryon/Photon Ratio 


T7 = (6.1±°;1) x lO" 10 



a Fit to the WMAP , CBI, ACBAR, 2dFGRS and Lyman a forest data 
b Assumes ionization fraction, x e = 1 

c I a = nd c /r s 
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can be incorporated in a single number: the optical depth due to reionization, r. The primordial 
fluctuations in this model are characterized by a spectral index. Despite its simplicity, it is an 
adequate fit not only to the WMAP temperature and polarization data but also to small scale 
CMB data, large scale structure data, and supernova data. This model is consistent with the 
baryon/photon ratio inferred from observations of D/H in distant quasars, the HST Key Project 
measurement of the Hubble constant, stellar ages and the amplitude of mass fluctuations inferred 
from clusters and from gravitational lensing. When we include large scale structure or Lyman a 
forest data in the analysis, the data suggest that we may need to add an additional parameter: 
dn s /d\xi k. Since the best fit models predict that the slope of the power spectrum is redder on small 
scales, this model predicts later formation times for dwarf galaxies. This modification to the power 
law ACDM model may resolve many of its problems on the galaxy scale. Table (10) lists the best 
fit parameters for this model. 

While there have been a host of papers on cosmological parameters, WMAP has brought this 
program to a new stage: WMAP 's more accurate determination of the angular power spectrum has 
significantly reduced parameter uncertainties, WMAP 's detection of TE fluctuations has confirmed 
the basic model and its detection of reionization signature has reduced the n s — r degeneracy. 
Most importantly, the rigorous propagation of errors and uncertainties in the WMAP data has 
strengthened the significance of the inferred parameter values. 

In this paper, we have also examined a number of more complicated models: non-flat universes, 
quintessence models, models with massive neutrinos, and models with tensor gravitational wave 
modes. By combining the WMAP data with finer scale CMB experiments and with other astro- 
nomical data sets (2dFGRS galaxy power spectrum and SNIa observations), we place significant 
new limits on these parameters. 

Cosmology is now in a similar stage in its intellectual development to particle physics three 
decades ago when particle physicists converged on the current standard model. The standard model 
of particle physics fits a wide range of data, but does not answer many fundamental questions: 
"what is the origin of mass? why is there more than one family?, etc." Similarly, the standard 
cosmological model has many deep open questions: "what is the dark energy? what is the dark 
matter? what is the physical model behind inflation (or something like inflation)?" Over the past 
three decades, precision tests have confirmed the standard model of particle physics and searched 
for distinctive signatures of the natural extension of the standard model: supersymmetry. Over 
the coming years, improving CMB, large scale structure, lensing, and supernova data will provide 
ever more rigorous tests of the cosmological standard model and search for new physics beyond the 
standard model. 
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Fig. 1. — This figure compares the best fit power law ACDM model to the WMAP temperature 
angular power spectrum. The gray dots are the unbinned data. 
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Fig. 2. — This figure compares the best fit power law ACDM model to the WMAP TE angular 
power spectrum. 
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Fig. 3. — This figure shows the likelihood function of the WMAP TT + TE data as a function 
of the basic parameters in the power law ACDM WMAP model. (Qbh 2 , Vt rn h?, h, A, n s and r.) 
The points are the binned marginalized likelihood from the Markov chain and the solid curve is an 
Edgeworth expansion of the Markov chains points. The marginalized likelihood function is nearly 
Gaussian for all of the parameters except for r. The dashed lines show the maximum likelihood 
values of the global six dimensional fit. Since the peak in the likelihood, xml is not the same as 
the expectation value of the likelihood function, < x >, the dashed line does not lie at the center 
of the projected likelihood. 
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Fig. 4. — This plot shows the contribution to 21n£ per multipole binned at Al = 15. The excess 
X 2 comes primarily from three regions, one around i ~ 120, one around I ~ 200 and the other 
around £ ~ 340. 
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Fig. 5. — Spectral Index Constraints. Left panel: the n s — r degeneracy in the WMAP data for a 
power-law ACDM model. The TE observations constrain the value of r and the shape of the Cf T 
spectrum constrain a combination of n s and r. Right panel: n s — f^/i 2 degeneracy. The shaded 
regions show the joint one and two sigma confidence regions. 
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Fig. 6. — (Left) This figure compares the best fit ACDM model of §3 based on WMAP data only to 
the 2dFGRS Power Spectrum(Percival et al. 2001). The bias parameter for the best fit Power Law 
ACDM model is 1.0 corresponding to a best fit value of (3 = 0.45. (Right) This figure compares the 
best fit Power Law ACDM model of §3 to the power spectrum at z = 3 inferred from the Lyman 
a forest data. The data points have been scaled downwards by 20%, which is consistent with the 
1 a calibration uncertainty (Croft et al. 2002). 
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Fig. 7. — (Left panel) This figure shows the fraction of mass in bound objects as a function of 
redshift. The black lines show the mass in collapsed objects with mass greater than M HRL (z), the 
effective Jeans mass in the absence of H2 cooling for our best fit PL ACDM model (thin lines are 
for the fit to WMAP only and thick lines are for the fit to all data sets). The heavy line uses the 
best fit parameters based on all data (which has a lower <rg) and the light line uses the best fit 
parameters based on fitting to the WMAP data only. The dashed lines show the mass in collapsed 
objects with masses greater than the Jeans mass assuming that the minimum mass is 1O 6 M . More 
objects form if the minimum mass is lower. (Right Panel) This figure shows the ionization fraction 
as a function of redshift. The solid line shows ionization fraction for the best fit PL ACDM model 
if we assume that H2 cooling is suppressed by photo-destruction of H2. This figure suggests that 
H2 cooling may be necessary for enough objects to form early enough to be consistent with the 
WMAP detection. The heavy line is for the best fit parameters for all data sets and the light line 
is for the best fit parameters for the WMAP only fit. The dashed lines assume that the objects 
with masses greater than 1O 6 M can form stars. The gray band shows the 68% likelihood region 
for z r based on the assumption of instantaneous complete reionization (Kogut et al. 2003). 
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Fig. 8. — This figure shows the marginalized likelihood for various cosmological parameters in the 
running spectral index model for our analysis of the combined WMAP , CBI, ACBAR, 2dFGRS 
and Lyman a data sets. The dashed lines show the maximum likelihood values of the global seven 
dimensional fit. 
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Fig. 9. — (Left)The shaded region in the figure shows the 1 — a contours for the amplitude of the 
power spectrum as a function of scale for the running spectral index model fit to all data sets. The 
dotted lines bracket the 2-a region for this model. The dashed line is the best fit power spectrum 
for the power law ACDM model. (Right)The shaded region in the figure shows the 1 -a contours 
for the amplitude of the amplitude of mass fluctuations, A 2 (k) = (k 3 /(27r 2 )P(k), as a function of 
scale for the running spectral index model fit to all data sets. The dotted lines bracket the 2-a 
region for this model. The dashed line is the best fit for the power law ACDM model. 
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Fig. 10. — (Left) This figure shows the fraction of the universe in bound objects with mass greater 
than M HRL (dashed), M j = 1O 6 M (solid) and M 3 '' (dotted) in a model with a running spectral 
index. The curves were computed for the la upper limit parameters for this model (see Figure 
9). These should be viewed as upper limits on the mass fraction in collapsed objects. (Right) 
This figure shows the ionization fraction as a function of redshift and is based on the assumptions 
described in §4.7. As in the figure on the left, we use the la upper limit estimate of the power 
spectrum so that we obtain "optimistic" estimates of the reionization fraction. In the context of a 
running spectral index fit to the data, the WMAP detection of reionization appears to require that 
i?2 cooling played an important role in early star formation. 
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Fig. 11. — Constraints on Dark Energy Properties. The upper left panel shows the marginalized 
maximum likelihood surface for the WMAPext data alone and for a combination of the WMAPext 
+ 2dFGRS data sets. The solid lines in the figure show the 68% and 95% confidence ranges for 
the fit C supernova data from Perlmutter et al. (1999). In the upper right panel, we multiply the 
supernova likelihood function by the WMAPext + 2dFGRS likelihood functions. The lower left 
panel shows the maximum likelihood surface for h and w for the WMAPext data alone and for the 
WMAPext + 2dFGRS data sets. The solid lines in the figures are the 68% and 95% confidence limits 
on Hq from the HST Key Project, where we add the systematic and statistical errors in quadrature. 
In the lower right panel, we multiply the likelihood function for the WMAPext + 2dFGRS data by 
the likelihood surface for the HST data to determine the joint likelihood surface. The dark areas 
in these plots are the 68% likelihood regions and the light areas are the 95% likelihood regions. 
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Fig. 12. — Constraints on Dark Energy Properties. The upper left panel shows the marginalized 
maximum likelihood surface for the WMAPext data alone and for a combination of the WMAPext 
+ 2dFGRS data sets. The solid lines in the figure show the 68% and 95% confidence ranges for 
supernova data from Riess et al. (2001). In the upper right panel, we multiply the supernova 
likelihood function by the WMAPext + 2dFGRS likelihood functions. The lower left panel shows 
the maximum likelihood surface for h and w for the WMAPext data alone and for the WMAPext 
+ 2dFGRS data sets. The solid lines in the figures are the 68% and 95% confidence limits on Hq 
from the HST Key Project, where we add the systematic and statistical errors in quadrature. In 
the lower right panel, we multiply the likelihood function for the WMAPext + 2dFGRS data by 
the likelihood surface for the HST data to determine the joint likelihood surface. The dark areas in 
these plots are the 68% likelihood regions and the light areas are the 95% likelihood regions. The 
calculations for this figure assumed a prior that w > —1. 
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Fig. 13. — Constraints on the geometry of the universe: VL m — Q\ plane. This figure shows the two 
dimensional likelihood surface for various combinations of data: (upper left) WMAP (upper right) 
WMAPext (lower left) WMAPext+ HST Key Project (supernova data (Riess et al. 1998, 2001) is 
shown but not used in the likelihood in this part of the panel; (lower right) WMAPext+ HST Key 
Project + supernova 
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Fig. 14. — This figure shows the marginalized cumulative probability of Q y h 2 based on a fit to 
the WMAPext+ 2dFGRS data sets (dashed) and the cumulative probability based on a fit to the 
WMAPext+ 2dFGRS+ Lyman a data sets (solid). The vertical lines are the 95% confidence upper 
limits for each case (0.21 and 0.23 eV). 
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Fig. 15. — This figure shows the cumulative likelihood of the combination of the WMAPext+ 
2dFGRS +Lyman a data sets as a function of r, the tensor/scalar ratio. The three lines show the 
likelihood for no priors, for models with \dn/dlnk\ < 0.005 and for models with n s < 1. 
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Fig. 16. — Angular correlation function of the best fit ACDM model, toy finite universe model, and 
WMAP data on large angular scales. The data points are computed from the template-cleaned V 
band WMAP using the KpO cut (Bennett et al. 2003c). 
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Fig. 17. — TE Power Spectrum. This figure compares the data to the predicted TE power spectrum 
in our toy finite universe model and the ACDM model. Both models assume that r = 0.17 and 
have identical cosmological parameters. This figure shows that the TE power spectrum contains 
additional information about the fluctuations at large angles. While the current data can not 
distinguish between these models, future observations could detect the distinctive TE signature of 
the model. 



